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(54) Signal Interpolation and decimation exploiting filter symmetry 



(57) Symmetry in a filter is used to reduce the com- 
plexity of an interpolator or a decimator and to simplify 
derivation of resulting discrete samples. In particular, an 
inverse relationship between weights applied to two 
samples is recognized and exploited. An inverse rela- 
tionship is recognized when a first weight is associated 
with a first of the samples and a second weight is asso- 
ciated with a second of the samples and a weight which 
is equivalent to the first weight is associated with the 
second sample and a weight which is equivalent to the 
second weight is associated with the first sample. The 
inverse relationship is exploited by forming two compos- 
ite weights of the first and second weights and weight- 
ing composite sample signals with the composite 
weights. A first of the composite weights has a value 
which is one-half of the sum of the values of the first and 
second weights. A second of the composite weights has 
a value which is one-half of the difference of the values 
of the first and second weights. The composite weights 
can be used repeatedly for each subsequent interpola- 
tion or decimation and are therefore calculated only 
once for processing many samples according to the 
same filter. The two composite samples have values 
which are, respectively, (i) the sum of the values of the 
first and second samples and (ii) the difference of the 
values of the first and second samples. Symmetry in 
weights applied to a single sample is recognized and 
exploited. For example, if at least two equivalent 
weights of a filter are applied to the same sample, one 
of the equivalent weights is applied to the sample once 
and the resulting weighted signal is used repeated 
rather than weighting the sample by each of the equiva- 
lent weights. Symmetry of weights applied to samples in 
producing a single resulting sample is recognized and 
exploited. For example, if equivalent weights are applied 



to two samples and the weighted samples are com- 
bined in producing a resulting sample, one of the equiv- 
alent weights is applied to the sum of the two samples. 
As a result, only one of the equivalent weights is applied 
to reduce the multiplication complexity of the interpola- 
tor or decimator. 
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Description 

FIELD OF THE INVENTION 

5 The present invention relates to signal processing in a computer system and, in particular, to a particularly efficient 
mechanism for interpolating and decimating discrete samples of an analog signal. 

BACKGROUND OF THE INVENTION 

10 Signal interpolation and signal decimation are forms of signal processing which require substantial resource in 
terms of computer resources and/or circuitry. Signal interpolation and decimation involve discrete samples of an analog 
signal generally taken at a specific frequency. Signal interpolation produces synthesized discrete signals of the analog 
signal at a greater frequency. For example, discrete samples of an analog audio signal taken at a frequency of 20MHz 
can be interpolated by a 1 -to-2 signal interpolator to produce synthesized discrete samples of the analog audio signal 

is at 40MHz. Signal decimation produces synthesized discrete signals of the analog signal at a lesser frequency. 

Signal interpolators and decimators are used to process discrete samples of various types of analog signals to pro- 
duce discrete samples of the analog signal at various frequencies. For example, signal interpolators and decimators 
can process discrete samples of an analog seismic signal recorded using a vibration source and one or more geo- 
phones in a conventional manner. Alternatively, signal interpolators and decimators can process discrete samples of an 

20 analog light signal, i.e., pixels of a color graphical image, recorded by a video camera. Furthermore, signal interpolators 
and decimators can process discrete samples of an analog audio recording. The discrete samples are typically taken 
from a source analog signal using an analog-to-digital converter, which converts a particular value of the analog signal 
at a particular time to a digital number which can be stored in the memory of a computer. 

Signal interpolation acid decimation typically requires significant amounts of processing resources due in part to 

25 the complexity of filters used in interpolation and decimation and due in part to the substantial number of discrete sam- 
ples processed. A filter is generally a number of weights which are applied to each of a number of discrete samples. 
The weights are generally referred to as a filter since the weights are applied to various collections of discrete samples. 
For example, a filter which has twenty-four weights is generally applied, first, to the first through twenty-fourth discrete 
samples; second, to the second through twenty-fifth discrete samples; third, to the third through twenty-sixth samples; 

30 and so on. In the case of a signal interpolator which has twenty-four weights, twenty-four multiplication operations and 
twenty-three addition operations are required to produce a synthesized discrete sample. Similarly, in the case of a sig- 
nal decimator which has twenty-four weights, twenty-four multiplication operations and twenty-three addition operations 
are required to process an original discrete sample of the analog signal. 

Signals which are interpolated or decimated typically include substantial numbers of discrete samples. For exam- 

35 pie, seismic data can involve scores of lines of seismic data, each line representing a path along the surface of the earth 
and each line including a thousand seismic traces or more. A seismic trace can tin turn include thousands of discrete 
samples of an analog seismic signal measured at a specific point on the surface of the earth. Therefore, a seismic sig- 
nal can easily include scores of millions of discrete samples. In another example, it is common today for graphical 
images to have one-thousand or more columns and one-thousand or more rows of picture elements, i.e., pixels, each 

40 of which is a discrete sample of a video signal. Such a video signal can be recorded using a video camera, an optical 
scanner, or can be generated by a computer to represent physical objects defined in part through physical manipulation 
of computer input devices by a user. Thus, it is common for a graphical image to include a million or more discrete sam- 
ples of a video image. In addition, motion video signals can include thousands of frames, each of which can include a 
million or more discrete samples Accordingly, interpolation or decimation of such seismic or video signals involves 

45 processing of millions of signals. Efficiency in a signal interpolator or decimator is therefore highly desirable to reduce 
the amount of time and resources required to process such signals. 

In addition, it is frequently desirable to interpolate or decimate signals very rapidly. For example, a compact disc 
player typically reads discrete samples of an analog audio signal at a rate of more than 40 million discrete samples per 
second. To enhance the sound quality of an analog audio signal reproduced from the discrete samples, additional dis- 

50 crete samples are interpolated from the discrete samples retrieved from the compact disc. Such an interpolator must 
generally process the discrete samples at a rate which is at least the rate at which the discrete samples are retrieved, 
i.e., at least 40 million discrete samples per second. 

Because of the significant processing resources required for such signal interpolation and decimation, a need per- 
sists in the industry for ever increasing efficiency in signal interpolators and decimators. 

55 

SUMMARY OF THE INVENTION 

In accordance with the present invention, symmetry in a filter is used to reduce the complexity of an interpolator or 
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a decimator and to simplify derivation of resulting discrete samples. In particular, an inverse relationship between 
weights applied to two samples is recognized and exploited. An inverse relationship is recognized when (i) a first weight 
is associated with a first of the samples and a second weight is associated with a second of the samples in the deriva- 
tion of a first resulting sample, and (ii) a weight which is equivalent to the first weight is associated with the second sam- 

5 pie and a weight which is equivalent to the second weight is associated with the first sample in a derivation of a second 
resulting sample. The inverse relationship is exploited by forming two composite weights of the first and second weights 
and weighting composite sample signals with the composite weights. A first of the composite weights has a value which 
is one-half of the sum of the values of the first and second weights. A second of the composite weights has a value 
which is one-half of the difference of the values of the first and second weights. The composite weights can be used 

10 repeatedly for each subsequent interpolation or decimation and are therefore calculated only once for processing many 
samples according to the same filter. The two composite samples have values which are, respectively, (i) the sum of the 
values of the first and second samples and (ii) the difference of the values of the first and second samples. Since only 
two composite weights are involved, two multipliers are required rather than four in which the first and second samples 
are weighted by the first and second weights, respectively, and by the second and first weights, respectively. The 

75 number of required adders is also reduced significantly as described more completely below. 

Further in accordance with the present invention, symmetry in weights applied to a single sample is recognized and 
exploited. For example, if two or more equivalent weights of a filter are applied to the same sample, one of the equiva- 
lent weights is applied to the sample once and the resulting weighted signal is used repeatedly rather than weighting 
the sample by each of the equivalent weights. 

20 Further in accordance with the present invention, symmetry of weights applied to samples in producing a single 
resulting sample is recognized and exploited. For example, if equivalent weights are applied to two samples and the 
weighted samples are combined in producing a resulting sample, one of the equivalent weights is applied to the sum of 
the two samples. As a result, only one of the equivalent weights is applied to reduce the multiplicative complexity of the 
interpolator or decimator. 

25 In general, the multiplicative complexity of an interpolator or decimator can be reduced by about 50% in accordance 
with the present invention if the weight filter is centrosymmetric. In addition, the additive complexity of an interpolator or 
decimator can be reduced asymptotically by 50%. Accordingly, interpolators and decimators according to the present 
invention represents a signal improvement in the efficiency of such interpolation and decimation. 

30 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram of a computer system which includes a signal processor in accordance with the present 
invention. 

Figures 2A-D illustrate a weight filter H as applied to source discrete samples X in interpolating samples Y. 
35 Figure 3 is a block diagram of a 1 -to-4 interpolator in accordance with the present invention. 
Figure 4 is a block diagram of a 1 -to-3 interpolator in accordance with the present invention. 
Figure 5 is a block diagram of a 3-to-1 decimator in accordance with the present invention. 
Figure 6 is a block diagram of a 1 -to-2 interpolator in accordance with the present invention in which additive com- 
plexity is reduced. 

40 Figure 7 is a block diagram of a 1 -to-5 interpolator in accordance with the present invention in which the weight filter 
is divided into two sub-filters. 

Figure 8 is a block diagram of a 1 -to-5 interpolator in accordance with the present invention in which additive com- 
plexity is reduced. 

45 DETAILED DESCRIPTION 

In accordance with the present invention, symmetry in a filter is used to reduce the complexity of an interpolator or 
a decimator and to simplify derivation of resulting discrete samples. 

so Hardware Components of the Imaoe Processing System 

To facilitate appreciation of the present invention, hardware components of one embodiment of a signal processing 
system in accordance with the present invention are briefly described. Computer system 100 (Figure 1) includes a proc- 
essor 1 02 and memory 1 04 which is coupled to processor 1 02 through a bus 1 06. Processor 1 02 fetches from memory 
55 104 computer instructions and executes the fetched computer instructions. Processor 102 also reads data from and 
writes data to memory 104 and sends data and control signals through bus 106 to one or more computer display 
devices 120 in accordance with fetched and executed computer instructions. Processor 102 can be, for example, any 
of the SPARC processors available from Sun Microsystems, Inc. of Mountain View, California. 
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Memory 104 can include any type of computer memory and can include, without limitation, randomly accessible 
memory (RAM), read-only memory (ROM), and storage devices which include storage media such as magnetic and/or 
optical disks. Memory 104 includes a signal processor 110, which is a computer process executing within processor 
102 from memory 104. A computer process is a collection of computer instructions and data which collectively define a 

5 task performed by computer system 1 00. As described more completely below, signal processor 1 1 0 (i) reads discrete 
source samples corresponding to a particular sampling frequency from a source buffer 112, (ii) derives from the dis- 
crete source samples new discrete resulting samples corresponding to a different sampling frequency, and (iii) stores 
those new discrete resulting samples in destination buffer 1 1 4. 

Each of computer display devices 120 can be any type of computer display device including without limitation a 

10 printer, a cathode ray tube (CRT), a light-emitting diode (LED) display, or a liquid crystal display (LCD). Each of compu- 
ter display devices 120 receives from processor 102 control signals and data and, in response to such control signals, 
displays the received data. Computer display devices 120, and the control thereof by processor 102, are conventional. 

Computer system 100 can be any of the SPARCstation workstation computers available from Sun Microsystems, 
Inc. of Mountain View, California. Sun, Sun Microsystems, and the Sun Logo are trademarks or registered trademarks 

75 of Sun Microsystems, Inc. in the United States and other countries. All SPARC trademarks are used under license and 
are trademarks of SPARC International, Inc. in the United States and other countries. Products bearing SPARC trade- 
marks are based upon an architecture developed by Sun Microsystems, Inc. 

Source buffer 1 12 and destination buffer 1 14 are stored in memory 104 and each contain discrete samples of an 
analog signal. For example, source buffer 1 1 2 can include discrete samples taken at a source sampling frequency from 

20 an analog seismic signal recorded using a vibration source and one or more geophones in a conventional manner. 
Alternatively, source buffer 112 can include discrete samples of an analog light signal, i.e., pixels of a color graphical 
image. Furthermore, source buffer 1 12 can include discrete samples taken at a source sampling frequency from an 
analog audio recording. The discrete samples of source buffer 1 1 2 are typically taken from a source analog signal using 
an analog-to-digital converter, which converts a particular value of the analog signal at a particular time to a digital 

25 number which can be stored in memory 104. The discrete samples are typically taken from the analog signal at a 
source sampling frequency such that the discrete samples represent values of the analog signal at regular time inter- 
vals. Destination buffer 114 generally includes discrete samples, which are derived by signal processor 1 10 from dis- 
crete source samples of source buffer 1 1 2 and which represent values of the analog signal corresponding to a different, 
destination sampling frequency. 

30 

1-to-N Interpolation 

Prior to describing specific efficiencies and benefits of the present invention, a brief discussion of interpolation is 
helpful. Figures 2A-D illustrate a 1 -to-4 interpolation. A weight filter H (Figure 2A) specifies relative weights for various 

35 discrete source samples of source buffer 1 1 2 (Figure 1 ). In this example, weight filter H (Figure 2A) includes 24 discrete 
weights HO-23. In Figure 2A, weight filter H is centered over destination sample Y0, which is to be interpolated from 
source samples X(N), X(N-1), X(N-2), X(N-3), X(N-4), and X(N-5). The value of destination sample Y0 is calculated by 
summing the values of source samples X(N), X(N-1), X(N-2), X(N-3), X(N-4), and X(N-5) as weighted by corresponding 
weights of weight filter H. As shown in Figure 2A, weights HO, H4, H8, H1 2, H1 6, and H20 are aligned with source sam- 

40 pies X(N), X(N-1), X(N-2), X(N-3), X(N-4), and X(N-5), respectively. Therefore, the value assigned to destination sample 
Y0 is determined by the following equation. 

Y0 = X(N) • HO + X(N-1) • H4 + X(N-2) • H8 + X(N-3) • H12 + X(N-4) • H16 + X(N-5) • H20 (1) 

45 Figure 2B shows the derivation of a second destination sample Y1 from source samples X(N), X(N-1), X(N-2), X(N- 
3), X(N-4), and X(N-5). In Figure 2B, weight filter H is moved to the left from the position of weight filter H in Figure 2A 
by an amount equal to the time interval between destination samples Y0 and Y1 . As shown in Figure 2B, weights H1 , 
H5, H9, H13, H17, and H21 are aligned with source samples X(N), X(N-1), X(N-2), X(N-3), X(N-4), and X(N-5), respec- 
tively. Therefore, the value assigned to destination sample Y1 is determined by the following equation. 

50 

Y1 = X(N) • H1 + X(N-1) • H5 + X(N-2) • H9 + X(N-3) • H1 3 + X(N-4) • H1 7 + X(N-5) • H21 (2) 

Figure 2C shows the derivation of a third destination sample Y2 from source samples X(N), X(N-1), X(N-2), X(N- 
3), X(N-4), and X(N-5). In Figure 2C, weight filter H is moved to the left from the position of weight filter H in Figure 2B 
55 by an amount equal to the time interval between destination samples Y1 and Y2. As shown in Figure 2C f weights H2, 
H6, H10, H14, H18, and H22 are aligned with source samples X(N), X(N-1), X(N-2), X(N-3), X(N-4), and X(N-5), 
respectively. Therefore, the value assigned to destination sample Y2 is determined by the following equation. 
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Y2 = X(N) • H2 + X(N-1) • H6 + X(N-2) • H10 + X(N-3) ♦ H14 + X(N-4) • H18 + X(N-5) • H22 (3) 

Figure 2D shows the derivation of a fourth destination sample Y3 from source samples X(N), X(N-1), X(N-2), X(N- 
3). X(N-4), and X(N-5). In Figure 2D, weight filter H is moved to the left from the position of weight filter H in Figure 2C 
by an amount equal to the time interval between destination samples Y2 and Y3. As shown in Figure 2D, weights H3, 
H7, H11, H15. H19, and H23 are aligned with source samples X(N), X(N-1), X(N-2), X(N-3), X(N-4). and X(N-5), 
respectively. Therefore, the value assigned to destination sample Y3 is determined by the following equation. 

Y3 = X(N) • H3 + X(N-1) • H7 + X(N-2) • H1 1 + X(N-3) • H15 + X(N-4) • H19 + X(N-5) • H23 (4) 

Weight filter H is typically symmetrical, i.e., weights HO, H1, H2, H3, H4, H5, H6, H7, H8, H9, H10, and H11 are 
equal to weights H23, H22, H21. H20, H19, H18, H17, H16. H15, H14, H13, and H12, respectively. In such a case, 
equations (1) through (4) are equivalent to equations (5) through (8), respectively. 

YO = X(N) • HO + X(N-1) • H4 + X(N-2) • H8 + X(N-3) • H1 1 + X(N-4) • H7 + X(N-5) • H3 (5) 

Y1 = X(N) ♦ H1 + X(N-1) • H5 + X(N-2) • H9 + X(N-3) • H10 + X(N-4) • H6 + X(N-5) • H2 (6) 

Y2 = X(N) • H2 + X(N-1) • H6 + X(N-2) • H10 + X(N-3) • H9 + X(N-4) • H5 + X(N-5) • H1 (7) 

Y3 = X(N) • H3 + X(N-1) • H7 + X(N-2) • H1 1 + X(N-3) • H8 + X(N-4) • H4 + X(N-5) • HO (8) 



Equations (5) through (8) is equivalent to the following equation (9). 



YO HO H4 H8 
Yl HI H5 H9 
Y2 = H2 H6 HIO 
Y3\ [H3 H7 Hll 



X(N) 

HI J H7 H3\ X(N-\) 
HIO H6 HI ,Y(iV-2) 
H9 H5 HI *(;V-3) 
H8 H4 H0\ ,Y(A'-4) 



Calculating interpolated values for Y0-Y3 in a straight forward, conventional manner according to equation (9) 
requires twenty-four multipliers and twenty adders. However, it is observed that weight HO is multiplied by sample X(N) 
and by sample X(N-5). In addition, weight H3 is also multiplied by samples X(N) and X(N-5). More specifically, HO • X(N) 
+ H3 • X(N-5) is a component of interpolated sample YO, and H3 • X(N) + HO • X(N-5) is a component of interpolated 
sample Y3. Together, these components require four (4) multipliers and two (2) adders. However, it is further observed 
that H0'X(N) + H3-X(N-5) is equivalent to (H0+H3)(X(N)+X(N-5)) + (H0-H3)(X(N)-X(N-5))/2. Similarly, H3-X(N) + 
H0-X(N-5)is equivalent to (H0+H3)(X(N)+X(N-5)) - (H0-H3)(X(N)-X(N-5))/2. The values (HO + H3)/2 and (HO - H3)/2 
recur in both components. Weights H0-H23 of weight filter H generally remain unchanged during the derivation of inter- 
polated samples of destination buffer 1 1 4 from all of the source samples of source buffer 112. Therefore, weights whose 
respective values are (HO + H3)/2 and (HO - H3)/2 can be precompiled and stored in registers within processor 102 or 
as data in memory 104 as composite weights cOO and dOO. 

Accordingly, components HO • X(N) + H3 • X(N-5) and H3 • X(N) + HO • X(N-5) are equivalent to (X(N) + X(N-5)) • 
cOO + (X(N) - X(N-5)) • dOO and (X(N) + X(N-5)) • cOO - (X(N) - X(N-5)) • dOO, respectively. Since (X(N) + X(N-5)) • 
cOO is included in both components, (X(N) + X(N-5)) • cOO is calculated once and used twice to form each of the two 
components. Similarly, (X(N) - X(N-5)) • dOO is included in both components and is calculated once and used twice to 
form each of the two components. Derivation of these components therefore requires two (2) multipliers and four (4) 
adders. 

At first glance, it appears that the overall complexity of an interpolator designed in accordance with equation (9) is 
not improved since four (4) multipliers in combination with two (2) adders is replaced with two (2) multipliers in combi- 
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nation with four (4) adders. However, two of the four (4) adders can be used to produce other components of interpo- 
lated signals produced in accordance with equation (9). For example, weight H1 is multiplied by sample X(N) and by 
sample X(N-5) and, weight H2 is also multiplied by samples X(N) and X(N-5). More specifically, H1 • X(N) + H2 • X(N-5) 
is a component of interpolated sample Y1, and H2»X(N) + H1 *X(N-5) is a component of interpolated sample Y2. 

5 H1^X(N) + H2-X(N-5) is equivalent to (H1+H2)(X(N)+X(N-5)) + (H1-H2)(X(N)-X(N-5))/2. Similarly, H2*X(N) + 
H1 • X(N-5) is equivalent to (H1 +H2)(X(N)+X(N-5)) - (H1 -H2)(X(N)-X(N-5))/2. The values (H1 + H2)/2 and (H1 - H2)/2 
recur in both components. Weights whose respective values are (H1+H2)/2 and (H1 - H2)/2 can be precompiled and 
stored in registers within processor 102 or as data in memory 104 as composite weights c10 and d10. Accordingly, 
components H1 • X(N) + H2 • X(N-5) and H1 • X(N) + H2 • X(N-5) are equivalent to (X(N) + X(N-5)) • c10 + (X(N) - X(N- 

w 5)) • d10 and (X(N) + X(N-5)) • c10 - (X(N) - X(N-5)) • d10, respectively. Since intermediate signals having respective 
values of (X(N) + X(N-5)) and (X(N) - X(N-5)) are produced in the production of interpolated samples Y0 and Y3 as 
described above, those same intermediate signals can also be used is producing interpolated signals Y1 and Y2. Spe- 
cifically, the production of the components of interpolated samples Y0-Y3 which involve source samples X(N) and X(N- 
5) collectively require four (4) multiplies and six (6) adders. 

15 In conventional interpolators designed in accordance with equation (9) in a straightforward manner require eight (8) 
multipliers and four (4) adders to produce components of interpolated samples Y0-Y3 which involve source samples 
X(N) and X(N-5). As described below with respect to Tables C and D, that interpolators and decimators in accordance 
with the present invention reduce the number of required multipliers to approximately 50% of the length of a weight filter 
H and reduce the number of required adders asymptotically to approximately 50% of the length of a weight filter H. 

20 Equation (9) above is therefore equivalent to the following equation. 

25 fro] [l 0 0 1 ] \c00 cOl c02 0 0 0 

Yl 0 i 1 0 clO ell c!2 0 0 0 

Y2 ~ 0 1 -1 0 0 0 0 d!2 dll dlO 

so Y3 1 G 0 -l] 1 0 0 0 ' d02 dOl dOO 



X{N)tX(N-5) 
X(N-\)+X{N-4) 
X(N-2)+X(N-3) 
X(N-2)-X(N-3) 
X(N-\)-X(N-4) 
X(N)-X(N-5) 



(10) 



35 

In equation (1 0), cij refers to a composite weight whose value is (H(i+4j) + H((3-i)+4j))/2 and dij refers to a compos- 
ite weight whose value is (H(i+4j) - H((3-i)44j))/2. For example, composite weight c12 has a value of (H9 + H10)/2, and 
composite weight d12 has a value of (H9 - Hl0)/2. 

It should be noted that the efficiencies described herein based upon the symmetry of weights HO and H23 and 

40 weights H3 and H20 with respect to source samples X(N) and X(N-5) can be achieved even if the remainder of weight 
filter H is not symmetrical. For example, the efficiencies described above can be achieved for any two symmetrical 
weights of weight filter H. Two weights, i.e., a first and a second weight, are symmetrical if the weights are equivalent 
and if, in calculating one interpolated sample, the first and second weights correspond to first and second source sam- 
ples, respectively, and, in calculating another interpolated sample, the first and second weights correspond to the sec- 

45 ond and first source samples, respectively. 

A four-to-one interpolator 300 (Figure 3) interpolates four interpolated samples Y0-Y3 from six source samples 
X(N) through X(N-5) according to equation (9) above. Source samples are received at terminal 302. At any given time, 
the signal received at terminal 302 represents source sample X(N). A series of delay units 304A-E preserve prior sig- 
nals received at terminal 302 to thereby record source samples X(N-1), X(N-2), X(N-3), X(N-4), and X(N-5) which are 

so produced on respective outputs of delay units 304A-E, respectively. Each of delay units 304A-E can be, for example, a 
shift register realized in processor 102 or in other circuitry or an item in a list stored in memory 104 (Figure 1). Negating 
circuitry 306C-E (Figure 3) receive from delay units 304C-E, respectively, source samples X(N-3), X(N-4), and X(N-5), 
respectively. Negating circuitry 306C-E produce signals which are the result of negating source samples X(N-3), X(N- 
4), and X(N-5), respectively. 

55 Interpolator 300 also includes adders 308A-F which produce signals representing the values (i) X(N) + X(N-5), (ii) 
X(N-1) + X(N-4), (iii) X(N-2) + X(N-3), (iv) X(N-2) - X(N-3), (v) X(N-1) - X(N-4), and (vi) X(N) - X(N-5), respectively. For 
example, adder 308A receives from terminal 302 source sample X(N) and from delay unit 302E source sample X(N-5) 
and produces a signal representing the sum of source samples X(N) and X(N-5). As a second example, adder 308F 
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receives from terminal 302 source sanrple X(N) and from negating circuitry 306E a signal representing the arithmetic 
negation of source sample X(N-5) and produces a signal representing the difference of source samples X(N) and X(N- 
5). 

Interpolator 300 further includes multipliers 310A-L Each of multipliers produces a signal which represents the 
arithmetic product of the signal produced by one of adders 308A-F and one of composite weights cOO, c01 , c02, c10, 
d 1 , c12, dOO, d01 , d02, d10. d1 1 , and d12 as shown. For example, multiplier 310A receives from adder 308A a signal 
representing the sum of source samples X(N) and X(N-5) and produces a signal which represents the arithmetic prod- 
uct of that sum and the value of composite weight cOO. Composite weight cOO is retrieved from memory 104 (Figure 1) 
or from a register of processor 1 02. As shown in Figure 3, the signal produced by adder 308A is used twice. Specifically, 
multiplier 310B also receives from adder 308A a signal representing the sum of source samples X(N) and X(N-5) and 
produces a signal which represents the arithmetic product of that sum and the value of composite weight c10. The val- 
ues represented by the respective signals produced by multipliers 310A-L are shown in Table A below. 



Table A 



Multiplier 


Value Represented By 
Produced Signals 


Source Signals From 


31 OA 


[X(N) + X(N-5)] • C00 


Adder 308A j 


310B 


[X(N) + X(N-5)] • c10 


Adder 308A 


310C 


[X(N-1) + X(N-4)] • c01 


Adder 308B 


310D 


[X(N-1) + X(N-4)J • c11 


Adder 308B 


310E 


[X(N-2) + X(N-3)] • C02 


Adder 308C 


310F 


[X(N-2) + X(N-3)J • C12 


Adder 308C 


31 0G 


(X(N-2)-X(N-3)] -d02 


Adder 308D 


310H 


[X(N-2)-X(N-3)] * d12 


Adder 308D 


3101 


[X(N-1)-X(N-4)] «d01 


Adder 308E 


310J 


[X(N-1)-X(N-4)J • d11 


Adder 308E 


310K 


[X(N)-X(N-5)] • dOO 


Adder 308F 


310L 


[X(N)-X(N-5)] -d10 


Adder 308F 



As shown in Table A and in Figure 3, the signal produced by each of adders 308A-F is used twice. Therefore, the 
number of adders required to produce interpolated samples Y0-Y3 is reduced from that required to produce interpo- 
lated signals using conventional techniques. In addition, no multiplication operations other than those performed by 
multipliers 310A-L are required to produce interpolated samples YO-3 as shown in Figure 3. Therefore, the number of 
multiplication operations is reduced by one-half, i.e., from twenty-four to twelve. The present invention therefore repre- 
sents a significant improvement over conventional interpolation techniques. 

As shown in Figure 3, interpolator 300 includes adders 31 2A and 31 2B which collectively produce a signal repre- 
senting the arithmetic sum of the signals produced by multipliers 310A, 310C, and 310E. Interpolator 300 also includes 
adders 31 2G and 312H which collectively produce a signal representing the arithmetic sum of the signals produced by 
multipliers 31 OH, 31 0J, and 310L. Adder 316A of interpolator 300 produces a signal representing the arithmetic sum of 
the signals produced by adders 312A and 312B and adders 312G and 312H. The signal produced by adder 316A is 
therefore interpolated sample Y0 according to equation (10) above. Similarly, adder 31 6D of interpolator 300 produces 
a signal representing the arithmetic difference of the signals produced by adders 31 2A and 31 2B and adders 31 2G and 
31 2H since adder 31 6D receives from negating circuitry 31 4B a signal representing the arithmetic negation of the signal 
produced by adders 312G and 312H. The signal produced by adder 316D is therefore interpolated sample Y3 accord- 
ing to equation (10) above. 

Interpolator 300 further includes adders 31 2C and 31 2D which collectively produce a signal representing the arith- 
metic sum of the signals produced by multipliers 310B, 310D, and 310E. In addition, interpolator 300 includes adders 
312E and 312F which collectively produce a signal representing the arithmetic sum of the signals produced by multipli- 
ers 31 0G, 3101, and 31 OK. Adder 316B of interpolator 300 produces a signal representing the arithmetic sum of the sig- 
nals produced by adders 312C and 312D and adders 312E and 312F. The signal produced by adder 316B is therefore 
interpolated sample Y1 according to equation (10) above. Similarly, adder 31 6C of interpolator 300 produces a signal 
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representing the arithmetic difference of the signals produced by adders 312C and 312D and adders 312E and 312F 
since adder 31 6C receives from negating circuitry 31 4A a signal representing the arithmetic negation of the signal pro- 
duced by adders 31 2E and 31 2F. The signal produced by adder 31 6C is therefore interpolated sample Y2 according to 
equation (10) above. 

In one embodiment, adders of interpolator 300, e.g., adders 308A-F, are computer instructions which, when exe- 
cuted by processor 1 02 (Figure 1 ) , perform addition operations in which operands are signals in the form of data stored 
in memory 104 of computer system 100. In this embodiment, multipliers of interpolator 300 (Figure 3), e.g., multipliers 
31 0A-L, are computer instructions which, when executed by processor 102 (Figure 1), perform multiplication operations 
in which operands are signals in the form of data stored in memory 104 of computer system 100. In an alternative 
embodiment, adders and multipliers of interpolator 300 (Figure 3) are voltage adders and voltage multipliers, respec- 
tively, and source samples, interpolated samples, and intermediate signals described above are voltages. In addition, 
interpolator 300 includes a number of instances of negation circuitry, e.g., negation circuitry 306E. In one embodiment, 
an adder, e.g, adder 308F, is capable of producing a signal which represents the arithmetic difference between respec- 
tive values of two received signals. In such an embodiment, negation circuitry such as negation circuitry 306E is obvi- 
ated and adders which receive signals produced by negation circuitry are configured to produce signals which 
represent arithmetic differences of the respective values or received signals. For example, negation circuitry 306E is 
omitted and adder 308F is configured to produce a signal representing the arithmetic difference between the respective 
values represented by signals received from terminal 302 and delay unit 304E in such an embodiment, in such an 
embodiment, negation circuitry does not add to the complexity of an interpolator or decimator. 

In the above illustrative example of the principles of the present invention, an even number of interpolated samples 
are derived from an even number of source samples. Further efficiencies can be realized when an odd number of inter- 
polated samples are derived or when interpolated samples are derived from an odd number of source samples. It is 
helpful to consider the following example. Suppose that three interpolated samples Y0-Y2 are derived from five source 
samples X(N) through X(N-4) from source buffer 112 (Figure 1). Equations which are analogous to equations (1) 
through (4) above can be written in matrix form as the following equation which describes the derivation of interpolated 
samples Y0-Y2. 











m ' 




Y0 




HO 


H3 


H6 H9 H12 






YJ 




HJ 


H4 


H7 H10 H13 


X{N-2) 


(U) 


Y2 




H2 


H5 


H8 Hli H14 














X(N-4) 





If the weight filter matrix is centrosymmetric, i.e., if weights HO, H1, H2, H3, H4, H5, and H6 are equivalent to 
weights H14, H13, H12, H1 1, H10, H9, and H8, respectively, equation (1 1) is equivalent to the following equation (12). 



Y0 




YJ 


— 


Y2 





X(N-2) 
,V(A'-3) 
X(N-4) 



(12) 



Efficiencies are realized in the manner described above based upon symmetry of weights HO, H2, H3, H5, H9, H 11 , 
H12, and H14. Further efficiency is realized from the symmetry of weights H1, H4, H10, and H13. 

Interpolated sample Y1 is derived by an interpolator in accordance with the value H1 • X(N) + H4 • X(N-1) + H7 
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• X(N-2) + H4 • X(N-3) + H1 • X(N-4) of equation (12). Efficiency is derived from multiple occurrences of weights H1 
and H4. Specifically, the value represented by interpolated sample is equivalent to H1 • [X(N) + X(N-4)] + H4 • [X(N- 
1) + X(N-3)] + H7 • X(N-2). Accordingly, by summing the respective values represented by source samples X(N) and 
X(N-4) with the respective values represented by source samples X(N-1) and X(N-3) prior to weighting the respective 
sums by weights H1 and H4. Accordingly, derivation of interpolated sample Y1 requires only three (3) multipliers and 
three (3) adders rather than five (5) multipliers and four (4) adders in a conventional interpolator. In general, summing 
source samples prior to weighting the sums with common weights as described reduces a calculation which ordinarily 
involves a number, n m , of multipliers and a number, n 3 =n m -1, of adders, the number of multipliers required are reduced 
to(n m +1)/2 and the number ofadders remains n 3 =n m -1. Interpolated sample Y1 which represents the value H1 • [X(N) 
+ X(N-4)] + H4 • [X(N-1) + X(N-3)] + H7 • X(N-2) is produced by a portion of an interpolator 400 (Figure 4) which imple- 
ments equation (12) in accordance the present invention. 

Interpolator 400 is described more completely below. Interpolator 400 is described briefly here and includes adders 
408A and 408B which produce respective signals representing the arithmetic sum of the values represented by source 
samples X(N) and X(N-4) and the arithmetic sum of the values represented by source samples X(N-1) and X(N-3), 
respectively. Interpolator 400 also includes multipliers 410A, 410C, and 410E. Multiplier 410A receives the signal pro- 
duced by adder 408A and weight H1 , which is scored in a register, in processor 1 02 or in memory 1 04 for example, and 
produces a signal which represents the arithmetic product of the signal produced by adder 408A and weight H1 . Simi- 
larly, multiplier 41 0C receives the signal produced by adder 408B and weight H4, which is also stored in a register in 
processor 1 02 or in memory 1 04 for example, and produces a signal which represents the arithmetic product of the sig- 
nal produced by adder 408B and weight H4. Multiplier 41 0E receives source sample X(N-2) from delay unit 404B and 
weight H7, which is also stored in a register in processor 102 or in memory 104 for example, and produces a signal 
which represents the arithmetic product of source sample X(N-2) and weight H7. The signals produced by multipliers 
41 OA, 41 0C, and 41 0E are then used to calculate interpolated samples Y0-Y2 in a manner described more completely 
below. 

Further efficiency is realized from the symmetry of weights H6 and H8 and from the position of weights H6 and H8 
in the same column of the weight filter matrix of equation (1 1). i.e., in that weights H6 and H8 are equivalent and applied 
to the same source sample in the derivation of separate interpolated samples. Specifically, from equation (12), compo- 
nents of the derivation of interpolated samples Y0 and Y2 both include H6 • X(N-2). A signal representing the value H6 

• X(N-2) is produced once by a multiplier 410F and used in deriving both interpolated sample Y0 and interpolated sam- 
ple Y2. The following equation, which is equivalent to equation (12) above, represents the efficiencies realized by inter- 
polator 400 in accordance with the present invention. 



Y0 




1 


0 


1 


( 


Yl 




0 


1 


0 




Y2 




1 


0 


-1 





Hi H4 H7 
0 0 0 



0 
0 



0 
0 



dOi dOO 



Y(A0+*(/V-4) 
X(N-\)+X(N-2) 

X{N-2) 
X(N-l)-X{N-3) 
X(N)-X{N-<\) 



(13) 



In equation (13), (i) the value represented by composite weight cOO is equivalent to (H0+H2)/2, (ii) the value repre- 
sented by composite weight c01 is equivalent to (H3+H5)/2, (iii) the value represented by composite weight dOO is 
equivalent to (H0-H2)/2, and (iv) the value represented by composite weight d01 is equivalent to (H3-H5)/2. Table B 
below shows the various connections and signals of interpolator 400 (Figure 4). 
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Table B 



Uuinpuneni 
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Y 
A 


V/KI\ 

X(N) 


r^a\a\i i inff AAA A 
Uclay Unil 4U4M 


Terminal A AO 

lerrninai 4Uc 


YfM 


Dplnv/ i mit AAAR 
Uolay unit *fU*f D 


Uelay Unil 4U4M 


YfhS 0\ 


Dplnv unit AC\AC\ 
L/ciay ui in tutu 


OpIpv/ 1 inrt AAAR 
Uclay Uilll HUHD 


A\y*-o) 


DpIav unit 4.04 D 


DpIsv unit ADAH 

UWCkJ UI III HUtv 




Monato ADfiP 
IMcyalo HUDw 


HpIow 1 mit AAAf* 
uciay unn qVrls 


-ANN -0) 


INcvjdlc HUDU 


PIqIqu 1 mit AHAn 
UWIdy Uilll ■H/HI^ 




AHHpr AnflA 


Tprminal AAO Dolov/ unit AAAH 
icmmidi fuc, ucidy uimtu*tu 




AHdar AAfiR 


DpIov/ 1 inrt AAA A Dplau i init AAA/"* 
UWIdy Uilll HuHn, Ucldy Uilll tWU 




AHHpr Anon 


DpIpv 1 rnrt AAAA Moaato AAfiP 
L/cldy Uilll HUta, licljclU? HUDvy 


YnM-1\ - Y(bl-0\ 
AUM 1 y " A^IM tj 


AHHpr AA.fi n 


lerrninai 4u^ f (Negate 4uou 


Y/M\ Y/M Vl\ 


N^i iliinlior A1AA 

iviuiTipiier 4 ium 


U1 AHHor AAflA 

n 1 , Muaer 4uom 


UI • rY/M\ 1 YJM >IYI 


iviuiiipiier 4 iub 


/«AA Arirtesr Af\QA 

cuu, Muaer 4UoM 


•an . rV/M\ 1 V/KI JIM 

COv • Ia(N) + A(N-4)J 


Multiplier 4 iuu 


n4, Muaer 4Uod 


M4 • |a(N-i ) + X(IM-2)J [ 


Mi i\i\r\\\ar A1An 

iviuriipuer 4 ivu 


rA1 AHHor AAfiR 

cui, Maaer wots 


rA1 • TY^M 1\ Y^M OM 
CUI •lA^IN-1; + A^IM-^JJ 


Mnltinlipr A10P 


H7 nplat/ 1 inrt AAAR 
rlf , Ucldy unn 4U4D 


H7 • Y/M-9\ 
n / • A^IM*<Sj 


Multiplier 41 OF 


H6, Delay unil 404B 


H6-X(N-2) 


Multiplier 41 OQ 


d01, Adder 408C 


d01«[X(N-1)-X(N-2)] 


Multiplier 41 OH 


dOO, Adder 408D 


d00«[X(N)-X(N-4)] 


Adder 41 2A 


Multiplier 41 OA, Multiplier 41 OC 


H1 • [X(N) + X(N-4)]+H4- [X(N-1) + X(N-2)] 


Adder 41 2B 


Adder 412A, Multiplier 41 OE 


H1 • [X(N) + X(N-4)]+H4 • [X(N-1 ) + X(N-2)]+H7 • X(N-2) = 
Yl 


Adder 41 2C 


Multiplier 410B, Multiplier 410D 


cOO • [X(N) + X(N-4)]+c01 ♦ [X(N-1) + X(N-2)] 


Adder 41 2D 


Adder 412C, Multiplier 41 OF 


COO • [X(N) + X(N-4)]+c01 • [X(N-1) + X(N-2)]+H6 • X(N-2) 


Adder 41 2E 


Multiplier 410G, Multiplier 410H 


d01 • [X(N-1) - X(N-2)]+dOO • p<(N) - X(N-4)] 


Negate 41 4E 


Adder 41 2E 


■d01 • [X(N-1) - X(N-2)]-d00 • [X(N) • X(N-4)] , 


Adder 41 6A 


Adder 41 2D, Adder 41 2E 


COO • [X(N) + X(N-4)]+c01 • [X(N-1 ) + X(N-2)J+H6 • X(N- 
2J-K101 • [X(N-1) - X(N-2)]-kI00 • [X(N) - X(N-4)] = YO 


Adder 41 6B 


Adder 412D, Negate 41 4E 


COO • [X(N) + X(N-4)]+c01 • [X(N-1) + X(N-2)]+H6 • X(N-2)- 
d01 • [X(N-1) - X(N-2)]-d00 • [X(N) - X(N-4)] = Y2 



Interpolator 400 (Figure 4) therefore uses eight (8) multipliers to produce interpolated samples Y0-Y2 whereas pro- 
ducing interpolated samples Y0-Y2 according to equation (11) above using conventional techniques requires fifteen 
(15) multipliers. In addition, interpolator 400 uses eleven (11) adders to produce interpolated samples Y0-Y2 whereas 
producing interpolated samples Y0-Y2 according to equation (1 1) above using conventional techniques requires twelve 
(12) adders. 

Decimators 

It is commonly known in digital network theory that an N-to-one decimator is the transpose of a one-to-N interpola- 
tor. Such is described, for example, in C.R. Crochiere and L.R. Rabiner, Multirate Digital Signal Processing . (Prentice 
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Hall, Englewood Cliffs, New Jersey) (1983). Accordingly, to design a particularly efficient N-toone decimator according 
to the present invention, a particularly efficient one-to-N interpolator is designed in accordance with the principles of the 
present invention and the transpose of the one-to-N interpolator is determined. 

For example, interpolator 400 (Figure 4) is a one-to-three interpolator in which weight filter H has a length of fifteen 
(15), i.e., includes fifteen (15) weights, and is centrosymmetric. A three-to-one decimator in which weight filter H has a 
length of fifteen (15) and is centrosymmetric is shown as decimator 500 (Figure 5). The design of decimator 500 is 
related to interpolator 400 (Figure 4) as follows. 

Each branch node of interpolator 400 in which a signal is duplicated and applied as an input signal to two compo- 
nents, such as adders or multipliers, is replaced with an adder whose input signals are supplied by the two components 
after transposition and whose output signal is the driver of the signal which is duplicated in the interpolator. For exam- 
ple, branch node 420D (Figure 4) of interpolator 400 duplicates the signal driven by delay unit 404D and applies the 
duplicated signal to negation circuitry 406D and to adder 408A. Transposition of branch node 420D results in adder 
520D (Figure 5) whose output signal is applied to delay unit 504D and whose input signals are driven by negation cir- 
cuitry 506D and adder 522A through branch 508A. Delay unit 504D and negation circuitry 506D correspond to delay 
unit 404D (Figure 4) and negation circuitry 406D, respectively. As described below, branch node 508A (Figure 5) is the 
result of transposition of adder 408A (Figure 4). 

Each adder of interpolator 400 is replaced with a branch node in decimator 500 (Figure 5) (i) whose input signal is 
driven by a component corresponding to components of interpolator 400 (Figure 4) which receives the output signal of 
the adder and (ii) whose output signals are applied to components corresponding to components of interpolator 400 
which supply the input signals of the adder. For example, the result of transposition of adder 408A is branch node 508A 
(Figure 5). Adder 408A (Figure 4) has input signals which are received from branch nodes 420A and 420D and an out- 
put signal which is supplied to branch node 422A. The results of transposition of branch nodes 420A, 420D, and 422A 
are adders 520A (Figure 5), 520D, and 522A, respectively. Accordingly, the result of transposition of adder 408A (Figure 

4) is branch node 508A (Figure 5) whose input signal is received from adder 522A and whose output signals are applied 
to adders 520A and 520D. 

Each of the remaining components of interpolator 400 (Figure 4) are transposed by swapping input signals and out- 
put signals. For example, multiplier 41 OH receives as an input signals the signal driven by adder 408D and multiplies 
the signal by composite weight dOO and supplies the resulting signal to adder 41 2E. The result of transposing multiplier 
41 OH is multiplier 51 OH (Figure 5) of decimator 500. Multiplier 51 OH (i) receives a signal from branch node 51 2E, which 
is the result of transposition of adder 41 2E (Figure 4), (ii) multiplies the received signal by composite weight dOO (Figure 

5) , and (iii) supplies the resulting signal to branch node 508D), which is the result of transposition of adder 408D (Figure 
4). 

The production of decimated sample components Y(N), Y(N-1), Y(N-2), Y(N-3), and Y(N-4) from source samples 
X0, X1 , and X2 by decimator 500 (Figure 5) is represented by the following equation (14). 







r(/V-i) 




Y(N-2) 


= 


Y{N-3) 




Y(N-4) 





1 0 0 0 1 

3 10 10 

0 0 I 0 0 

0 10-1 0 

1 0 0 0 -1 



cOO HI 0 

cOI H4 0 

H6 H7 0 

0 0 dOl 
0 0 dOO 



X0+X2 

XI 
X0-X2 



(14) 



Decimator 500 (Figure 5) is the result of transposition of interpolator 400 (Figure 4). The efficiencies realized in 
interpolator 400 as described above are similarly realized in decimator 500 (Figure 5). Analogously, any of the efficien- 
cies described herein with respect to an interpolator can also be realized in a decimator which is related to an interpo- 
lator which is designed in accordance with the principles of the present invention. For example, composite weights cOO 
and c01 are applied to signals representing the sum of the values represented by source samples X0 and X2 and com- 
posite weights dOO and d01 are applied to signals representing the difference of the values represented by source sam- 
ples X0 and X2 in a directly analogous manner to that described above with respect to interpolator 300 (Figure 3). In 
addition, weights H1 and H4 are applied to source sample X1 once each and used twice each in the production of dec- 
imated sample components Y(N), Y(N-1), Y(N3), and Y(N-4) in a directly analogous manner to that described above 
with respect to interpolator 400 (Figure 4). Also in a directly analogous manner to that described above with respect to 
interpolator 400, weight H6 is applied to a signal representing the sum of the values represented by source samples X0 



12 



EP0818740A2 



and X2 within decimator 500 (Figure 5). 
Passing Qgrnplexity 

A useful measurement of the efficiency of an interpolator or decimator is the processing complexity. With respect 
to an interpolator, processing complexity represents a number of processing components, e.g., multipliers and adders, 
for each input. For example, the processing complexity of interpolator 400 (Figure 4) represents the number of process- 
ing components required to produce interpolated samples Y0, Y1 , and Y2 before source samples X(N) through X(N-4) 
are shifted through delay units 404A-D to read another source sample at terminal 402. With respect to a decimator, 
processing complexity represents a number of processing components for each output. For example, the processing 
complexity of decimator 500 (Figure 5) represents the number of processing components required to produce a deci- 
mated sample at terminal 502 from three source samples received at terminals 528A-C. 

Table C below shows the processing complexity of a one-to-N interpolator designed in accordance with the princi- 
ples of the present invention and which involves a weight filter H which has a number, L, of weights and which is cen- 
trosymmetric. In Table C below, M = L7N and N, L, and M are integers. 



Table C 



M 


N 


Multipliers 


\ Adders 


Delay units 


even 


even 


L/2 


L/2+M 


M-1 


even 


odd 


L/2 


U2+M-1 


M-1 


odd 


even 


L/2 


LV2+M-1 


M-1 


odd 


odd 


(L+D/2 


(L+1)/2+M-2 


M-1 



By comparison, generating interpolated signals Y0-Y2 according to equation (1 1) above requires generally L mul- 
tipliers and L-N adders. 

Table D below shows the processing complexity of an N-to-one decimator designed in accordance with the princi- 
ples of the present invention and which involves a weight filter H which has a number, L, of weights and which is cen- 
trosymmetric. In Table D below, M ° UN and N, L, and M are integers. 



Table D 



M 


N 


Multipliers 


Adders 


Delay units 


even 


even 


172 


L/2+M+N-1 


M-1 


even 


odd 


L/2 


L/2+M+N-2 


M-1 


odd 


even 


L/2 


L/2+M+N-2 


M-1 


odd 


odd 


(L+1)/2 


(L+1)/2+M+N-3 


M-1 



By comparison, a decimator which is the transpose of a conventional interpolator designed according to equation 
(1 1) above generally L multipliers and L-1 adders. 

Fgrther Improvement in Active Complexity 

Additive complexity, i.e., the number of adders of an interpolator or decimator in accordance with the present inven- 
tion, can be further reduced in some circumstances. TTie following is an illustrative example. 

The following equation represents the interpolation of two interpolated samples using a weight filter H whose length 
is twenty-four (24). 
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YO 




HO 


H2 


H4 


H6 


H8 


mo 


H12 


H14 


H16 


H18 


Yl 




P 1 


H3 


H5 


H7 


H9 


Hll 


H13 


HIS 


H17 


H19 



m 
w 

X{N-2) 
X(M-3) 
X(N-4) 
X[N-5) 
X(N-6) 
X{N-1) 
X{N-%) 
XlN-9) 
*(/V-lO) 
X(N-\\) 



The additive complexity of an interpolator designed according to equation (15) is twenty-four (24) as can be seen 
in Table C above. In other words, for each interpolation to produce two interpolated samples YO and Y1, twenty-four 
adders are required. 

The weight matrix of equation (15) has twelve columns and two rows. By recognizing that source sample X(N-1) is 
source sample X(N) after one time delay unit in the sampling rate of source buffer 1 1 2 (Figure 1 ), equation (1 5) can be 
rewritten as the following equation (16). 



YO 




Yl 





1 : x ' 1 0 0 
0 0 1 :{■* 



HO H4 H8 H12 H16 H20 

H2 H6 HIO HN H18 H22 

HI H5 H9 HI 3 HJ7 H21 

H3 H7 Hll H15 H19 H23 



X(N) 
X(N-2) 

X(N-S) 
IA'(/V-10) 



(16) 



In equation (16), z-f 1 represents a single time delay in the sampling rate of source buffer 1 12. In other words, inter- 
polated sample YO is produced by (i) weighting each of source samples X(N), X(N-2), X(N-4), X(N-6), X(N-8), and X(N- 
1 0) with a respective weight of the first row of the weight matrix of equation (1 6), (ii) weighting each of the same source 
samples with a respective weight of the second row of the weight matrix of equation (16) after a single delay, and (iii) 
combining the source samples as weighted to produce interpolated sample YO. 

Interpolator 600 (Figure 6) realizes the interpolation of equation (16). Interpolator 600 includes circuitry 602 which 
implements a one-to-four interpolator which is generally analogous to interpolator 300 (Figure 3) as described above. 
Interpolator 600 (Figure 6) includes delay units 604A-E each of which delays the source signal received at terminal 610 
by an amount equal to twice the sampling interval of source buffer 1 12 (Figure 1). Therefore, whereas interpolator 300 
(Figure 3) processes source samples X(N), X(N-1), X(N-2), X(N-3), X(N-4), and X(N-5). interpolator 600 (Figure 6) 
processes source samples X(N), X(N-2), X(N-4), X(N-6), X(N-8), and X(N-10). 

In a manner analogous to that described above with respect to interpolator 300 (Figure 3), circuitry 602 (Figure 6) 
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of interpolator 600 produces four output signals on each of terminals 612A-D and each of which corresponds to a 
respective row of the weight matrix of equation (16). To produce interpolated sample Y0 in the manner described above, 
the produced signal on terminal 61 2B is delayed by delay unit 606A for an amount of time equal to a single sampling 
interval of source buffer 112 (Figure 1). The delayed signal is added to the signal produced on terminal 61 2A by adder 
608A. The resulting signal is interpolated sample Y0. Interpolator 600 produces interpolated sample Y1 is a directly 
analogous manner using terminals 612C-D, delay unit 608B, and adder 608B. 

The additive complexity of interpolator 600 is twenty (20): eighteen (18), which is the additive complexity of circuitry 
602 according to Table C above, plus two (2), as a result of adders 608A-B. Thus, four (4) fewer adders a required by 
interpolator 600 than to implement an interpolator according to equation (15) above. 

Weight Filters of Arbitrary length 

In the illustrative examples described above, the number of weights of the weight filter H is an integer multiple of 
the number of interpolated samples in the context of interpolators or an integer multiple of the number of source sam- 
ples in the context of decimators. However, the principles of the present invention are equally applicable to interpolators 
and decimators in which the number of weights of the weight filter H is not an integer multiple of the number of interpo- 
lated samples of an interpolator or an integer multiple of the number of source samples of a decimator. The following 
example is illustrative. 

The following equation (17) corresponds to interpolation of seven (7) interpolated samples using a weight filter 
which includes thirty-three (33) weights. 
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(17) 



In this illustrative example, weights H0-H32 are centrosymmetric. Accordingly, equation (1 7) is equivalent to the fol- 
lowing equation (18). 
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Equation (18) can be divided into the following two separate equations (19) and (20). 
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(19) 
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The weights of the weight filter matrices of equations (1 9) and (20) are centrosy mmetric and can therefore each be 
independently implemented in the manner described above to achieve the efficiencies described above with respect to 
interpolators 300 (Figure 3) and 400 (Figure 4) for example. 

Interpolator 700 (Rgure 7) interpolates seven (7) interpolated samples using thirty-three (33) weights in accord- 
ance with equations (1 7)-(20). Interpolator 700 includes delay units 702A-D to form source samples X(N), X(N-1 ), X(N- 
2), X(N-3), and X(N-4) from source sample buffer 112 (Figure 1) in the manner described above. Source samples X(N), 
X(N-1), X(N-2), X(N-3), and X(N-4) are applied to circuitry 704 which produces interpolated samples Y0, Y1, Y2, Y3 P 
and Y4 in accordance with equation (19) in a manner which is analogous to that described above with respect to inter- 
polators 300 (Figure 3) and 400 (Figure 4). Source samples X(N), X(N-1), X(N-2), and X(N-3) are applied to circuitry 
706 which produces interpolated samples Y5 and Y6 in accordance with equation (20) in a manner which is analogous 
to that described above with respect to interpolators 300 (Figure 3) and 400 (Figure 4). 

Thus, efficiencies described above in interpolators and decimators which involve a number of weights which is an 
integer multiple of the number of output sample per interpolation or input samples per decimation can be realized inter- 
polators and decimators in which the number of weights is not an integer multiple of the number of output sample per 
interpolation or input samples per decimation. 

Reducing Additive Complexity in Conjunction with Weight Filters of Arbitrary Length 

Interpolator 700 (Figure 7) represents a substantial improvement in efficiency over conventional interpolators. It is 
appreciated, however, that circuitry 704 and circuitry 706 share input signals, e.g., source samples X(N), X(N-1), X(N- 
2), X(N-3), but share no signals subsequently derived from the input signals. Since source samples X(N) and X(N-2) 
are equivalent to source samples X(N-1) and X(N-3), respectively, after one time delay in the sampling rate of source 
buffer 112 (Figure 1), equation (18) above is equivalent to equation (21) below, in which z-f 1 represents a single time 
delay in the sampling rate of source buffer 112. 
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(21) 



Derivation of interpolated samples Y0-Y4 includes intermediate signals corresponding to values X(N-1)+X(N-3) 
and X(N-1 )-X(N-3) for generally the reasons described above with respect to equations (9) and (10). According to equa- 
tion (21), the same intermediate signals are used to derive interpolated signals Y5 and Y6. Specifically, weights H5 and 
H6 which are applied to source samples X(N) and X(N-3) according to equations (18) and (20) are applied to source 
samples X(N-1) and X(N-3) in equation (21). In addition, weights H12 and H13 which are applied to source samples 
X(N-1) and X(N-2) according to equations (18) and (20) are applied to source samples X(N-1) and X(N-3) in equation 
(21). As a result, weights H5, H6, H12, and H13 are applied to source samples X(N-1) and X(N-3) in equation (21) as 
are weights H7-H11. Accordingly, symmetry in the remainder the matrix representation of weight filter H is further 
exploited by aligning the symmetrical application of weights H5, H6, H12, and H13 with the symmetrical application of 
weights H7-H1 1 of weight filter H. 

Equation (21) is equivalent to equation (22) below in which interpolated signals Y5* and Y6' are equivalent to inter- 
polated signals Y5 and Y6, respectively, corresponding to a set of source samples X(N), X(N-1), X(N-2), X(N-3), and 
X(N-4) before one time delay in the sampling rate of source buffer 112 (Figure 1). 
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Interpolator 800 (Figure 8) implements the derivation of interpolated signals Y5\ Y6\ YO, Y1, Y2, Y3, and Y4 
according to equation (22). Interpolator 800 includes delay units 802A-D which are directly analogous to delay units 
702A-D (Figure 7) and which separate source signals X(N), X(N-1), X(N-2), X(N-3), and X(N-4) (Figure 8) in a directly 
analogous manner to that described above with respect to delay units 302 A-E (Figure 3). Interpolator 800 (Figure 8) 
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includes circuitry 804 which is directly analogous to circuitry 704 (Figure 7) and which produces interpolated samples 
Y0, Y1 , Y2, Y3, and Y4 (Figure 8) in accordance with equation (1 9) in a analogous manner to that described above with 
respect to interpolators 300 (Figure 3) and 400 (Figure 4). 

Intermediate signals X(N-1)+X(N-3) and X(N-1)-X(N-3) of circuitry 804 are used to derive interpolated signals Y5* 
and Y6' according to equation (22) above. Table E below shows the various connections and signals of interpolator 800 
(Figure 8) which operate to derive interpolated signals Y5' and Y6' In interpolator 800, and therefore Table E below, 
composite weight aOO is equivalent to (H5+H13)/2, composite weight a10 is equivalent to (H6+H12)/2, composite 
weight bOO is equivalent to (H5-H13)/2, and composite weight b10 is equivalent to (H6-H12)/2. 



Table E 



Component 


Signals Received From 


Value Represented by Output Signal 


Multiplier 806A 


a10, Circuitry 804 


a10-[X(N-1) + X(N-3)] 


Multiplier 806B 


aOO, Circuitry 804 


a00-[X(N-1) + X(N-3)] 


Multiplier 806C 


b10, Circuitry 804 


b10-[X(N-1)-X(N-3)] 


Mi il+irxliar* QHCPl 

Multiplier oudu 


duu. uircuitry ou4 


DOU • Ia(IM-I) • a(N-3)J 


Negate 808C 


Multiplier 806C 


-b10«[X(N-1)-X(N-3)] 


Negate 808D 


Multiplier 806D 


-b00«[X(N-1)-X(N-3)] 


Adder 81 OA 


Multiplier 806A, Multiplier 806C 


a10 • [X(N-1) + X(N-3)] + b1 0 • [X(N-1) - X(N-3)] 


Adder 810B 


Multiplier 806B, Multiplier 806D 


a00 • [X(N-1) + X(N-3)] + bOO • [X(N-1) - X(N-3)] 


Adder 81 0C 


Multiplier 806A, Negate 808C 


a10 • [X(N-1) + X(N-3)] - b10 • [X(N-1) - X(N-3)] 


Adder 81 0D 


Multiplier 806B, Negate 808D 


aOO • [X(N-1) + X(N-3)] + bOO • [X(N-1) - X(N-3)] 


Delay unit 81 2C 


Adder 810C 


z,- 1 • [a10 • [X(N-1) + X(N-3)] - b10 • [X(N-1) - X(N-3)]] 


Delay unit 81 2D 


Adder 810D 


z, -1 • [aOO • [X(N-1) + X(N-3)] + bOO • [X(N-1) - X(N-3)]] 


Adder 814A 


Adder 810A, Delay unit 812C 


a10 • [X(N-1) + X(N-3)] + b10 • [X(N-1) - X(N-3)] + z^ 
1 • [a10 • [X(N-1) + X(N-3)] - b1 0 • [X(N-1) - X(N-3)]] =Y6' 


Adder 81 4B 


Adder 81 0B, Delay unit 81 2D 


aOO • [X(N-1) + X(N-3)] + bOO • [X(N-1) - X(N-3)] + z{ 
1 • [aOO • [X(N-1) + X(N-3)] + bOO • [X(N-1) - X(N-3)D = Y5' 



In producing interpolated signals Y5' and Y6', interpolator 800 (Figure 8) includes four (4) multipliers and six (6) 
adders. By comparison, circuitry 706 (Figure 7) of interpolator 700 includes four (4) multipliers and nine (9) adders to 
produce interpolated signals Y5 and Y6 as shown in Table D. Accordingly, interpolator 800 (Figure 8) produces interpo- 
lated signals Y5 and Y6 with three (3) fewer adders by constructing interpolator 800 according to equation (22) in which 
symmetrical application of weights H5 f H6, H12, and H13 is aligned with symmetrical application of weights H7-H11 
such that each of weights H5-H13 is applied to both source samples X(N-1) and X(N-3). 

The additive complexity of interpolators in which a weight function H has an arbitrary length is shown below in Table 
F. In Table F, weight function H has a length L which is equal to MN+K where N is the number of interpolated signals 
produced by the interpolator, and W equals M plus the integer portion of (L+1)/2. L, M, N, K, and W are integers. 
Approach A refers to the general approach described above with respect to interpolator 700 (Figure 7) and equations 
(18) through (20), and Approach B refers to the general approach described above with respect to interpolator 800 (Fig- 
ure 8) and equations (21) and (22). 
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Table F 
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Additive Complexity 
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20 In general, Approach A and Approach 6 require the same number of multipliers, but Approach B requires additional 
delay units. Specifically, Approach B requires N-K additional delay units when M is even and K additional delay units 
when M is odd. However, since delay units are relatively inexpensive, it is preferred that the one of Approaches A and 
B which has the least additive complexity is used. 

25 Nth-Band Filters 

Nth-band filters are commonly used for 1-to-N interpolators and for N-to-1 decimators and are sometimes referred 
to as Nyquist filters. An Nth-band filter has a center weight, i.e., H c , and therefore has an odd length. In addition, an 
Nth-band filter includes weights H c+Ni , whose values are zero for each positive or negative integer, i. The signal sample 
30 interpolation and decimation mechanisms described above are directly applicable to increasing efficiency of interpola- 
tors and decimators which use Nth-band filters. 

The following example is illustrative. Consider a 4th-band interpolator with a weight filter H of 19 weights. The fol- 
lowing equation (23) represents the derivation of interpolated signals Y0-Y3. 
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45 



The center weight is weight H9, and weight filter H is centrosymmetric. Accordingly, equation (23) is equivalent to 
equation (24) below. 
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In equation (24), weight filter H has a length of L which is NM+K where N=4, M=4, and K=3. An interpolator can 
therefore be designed in accordance with equation (24) in the manner described above with respect to interpolators 700 
(Figure 7) and 800 (Figure 8). In addition, weights H1 and H5 are equal to zero since N=4 and H 9±Nj =0 for each integer, 
i. Accordingly, equation (24) is equivalent to equation (25) below. 
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It should be noted that the first and third rows of weight matrix H are symmetrical with respect to one another and 
therefore readily lend themselves to the improvements in efficiency described above with respect to interpolators 300 
(Figure 3) and 400 (Figure 4). In addition, the second and fourth rows of weight matrix H are each individually symmet- 
rical and therefore readily lend themselves to improvements in efficiency described above with respect to interpolator 
400 (Figure 4). Furthermore, weight H8 is applied to source sample X(N-2) in deriving both interpolated samples Y0 
and Y2. Accordingly, analogous improvements in efficiency described above with respect to interpolator 400 (Figure 4) 
are also applicable to an interpolator designed in accordance with equation (25). 

The above description is illustrative only and is not limiting. The present invention is limited only by the claims which 
follow. 

Claims 

1 . A method for deriving at least one resulting samples of a resulting signal from at least two source samples of a 
source signal, the method comprising: 

associating each of at least five weights with one of the at least two source samples; 
determining that the respective weight associated with a first of the source samples in the derivation of a first 
of the at least two resulting samples is substantially equivalent to the respective weight associated with a sec- 
ond of the source samples in the derivation of a second of the at least two resulting samples; 
determining that the respective weight associated with the second source sample in the derivation of the first 
resulting sample is substantially equivalent to the respective weight associated with the first source sample in 
the derivation of the second resulting sample; 

forming first and second composite weights of the respective weights associated with the first and second 

source samples in the derivation of the first resulting sample; 

storing the first and second composite weights in the memory of a computer; 

forming a sum signal from the first and second source samples; 

forming a difference signal from the first and second source samples; 

weighting the sum signal by the first composite weight; and 
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weighting the difference signal by the second composite weight. 

2. The method of Claim 1 further comprising: 

5 determining that a selected of the weights is associated with a selected one of the source samples in the der- 

ivation of a third of the resulting samples; 

determining that another of the weights, which is equivalent to the selected weight, is associated with the 
selected source sample in the derivation of a fourth of the resulting samples; 
weighting the selected source sample with the selected weight to form a selected weighted signal; and 
10 including the selected weighted signal in deriving both the third and the fourth resulting samples. 

3. The method of Claim 2 wherein the third resulting sample is the first resulting sample. 

4. The method of Claim 1 further comprising: 

15 

determining that a selected one of the weights is associated with a third one of the source samples in the der- 
ivation of a third one of the resulting samples; 

determining that another of the weights, which is equivalent to the selected weight, is associated with a fourth 
one of the source samples in the derivation the third resulting sample; 
20 forming a sum signal from the third and fourth source samples; 

weighting the sum signal with the selected weight to form a weighted sum signal; and 
including the weighted sum signal in deriving the third resulting sample. 

5. The method of Claim 4 wherein the third source sample is the first source sample. 

25 

6. The method of Claim 4 wherein the fourth source sample is the second source sample. 

7. A computer program product comprising a computer usable medium having computable readable code embodied 
therein for deriving at least one resulting samples of a resulting signal from at least two source samples of a source 

30 signal, the computer readable code comprising: 

at least five weights each of which are associated with one of the at least two source samples, 

wherein the respective weight associated with a first of the source samples in the derivation of a first of the 
35 at least two resulting samples is substantially equivalent to the respective weight associated with a second 

of the source samples in the derivation of a second of the at least two resulting samples, and 
further wherein the respective weight associated with the second source sample in the derivation of the 
first resulting sample is substantially equivalent to the respective weight associated with the first source 
sample in the derivation of the second resulting sample; 

40 

a composite weight builder operatively coupled to the at least five weights and confirmed to form first and sec- 
ond composite weights of the respective weights associated with the first and second source samples in the 
derivation of the first resulting sample; 

an sum adder configured to form a sum signal from the first and second source samples; 
45 a difference adder configured to form a difference signal from the first and second source samples; 

a sum multiplier operatively coupled to the sum adder and configured to weight the sum signal by the first com- 
posite weight; and 

a difference multiplier operatively coupled to the difference adder and the composite weight builder and config- 
ured to weight the difference signal by the second composite weight. 

50 

8. The computer program product of Claim 7 wherein the at least five weights include a selected weight associated 
with a selected one of the source samples in the derivation of a third of the resulting samples and another weight 
equivalent to the selected weight and associated with the selected source sample in the derivation of a fourth of the 
resulting samples; 

55 

further wherein the computer readable code further comprises: 

a selected multiplier configured to weight the selected source sample with the selected weight to form a 
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selected weighted signal; 

a third resulting sample builder operatively coupled to the selected multiplier and configured to produce the 
third resulting sample using the selected weighted signal; and 

a fourth resulting sample builder operatively coupled to the selected multiplier and configured to produce 
5 the fourth resulting sample using the selected weighted signal. 

9. The computer program product of Claim 8 wherein the third resulting sample is the first resulting sample. 

10. The computer program product of Claim 7 wherein the at least five weights include a selected weight associated 
w with a third one of the source samples in the derivation of a third one of the resulting samples and another of the 

weights equivalent to the selected weight and associated with a fourth one of the source samples in the derivation 
the third resulting sample; 

further wherein the computer readable code further comprises: 

15 

a selected adder configured to form a sum signal from the third and fourth source samples; 

a selected multiplier operatively coupled to the selected adder and configured to weight the sum signal 

with the selected weight to form a weighted sum signal; and 

a third resulting sample builder operatively coupled to the selected multiplier and configured to produce the 
20 third resulting sample using the weighted sum signal. 

1 1 . The computer program product of Claim 1 0 wherein the third source sample is the first source sample. 

12. The computer program product of Claim 10 wherein the fourth source sample is the second source sample. 

25 

13. A signal processor comprising: 

at least five weights each of which are associated with one of at least two source samples, 

30 wherein the respective weight associated with a first of the source samples in the derivation of a first of at 

least two resulting samples is substantially equivalent to the respective weight associated with a second 
of the source samples in the derivation of a second of the at least two resulting samples, and 
further wherein the respective weight associated with the second source sample in the derivation of the 
first resulting sample is substantially equivalent to the respective weight associated with the first source 

35 sample in the derivation of the second resulting sample; 

a composite weight builder operatively coupled to the at least five weights and configured to form first and sec- 
ond composite weights of the respective weights associated with the first and second source samples in the 
derivation of the first resulting sample; 
40 an sum adder configured to form a sum signal from the first and second source samples; 

a difference adder configured to form a difference signal from the first and second source samples; 
a sum multiplier operatively coupled to the sum adder and configured to weight the sum signal by the first com- 
posite weight; and 

a difference multiplier operatively coupled to the difference adder and the composite weight builder and config- 
45 ured to weight the difference signal by the second composite weight. 

14. The signal processor of Claim 13 wherein the at least five weights include a selected weight associated with a 
selected one of the source samples in the derivation of a third of the resulting samples and another weight equiv- 
alent to the selected weight and associated with the selected source sample in the derivation of a fourth of the 

so resulting samples; 

further wherein the signal processor further comprises: 

a selected multiplier configured to weight the selected source sample with the selected weight to form a 
55 selected weighted signal; 

a third resulting sample builder operatively coupled to the selected multiplier and configured to produce the 
third resulting sample using the selected weighted signal; and 

a fourth resulting sample builder operatively coupled to the selected multiplier and configured to produce 
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the fourth resulting sample using the selected weighted signal. 

15. The signal processor of Claim 14 wherein the third resulting sample is the first resulting sample. 

5 1 6. The signal processor of Claim 1 3 wherein the at least five weights include a selected weight associated with a third 
one of the source samples in the derivation of a third one of the resulting samples and another of the weights equiv- 
alent to the selected weight and associated with a fourth one of the source samples in the derivation the third result- 
ing sample; 

10 further wherein the signal processor further comprises: 

a selected adder configured to form a sum signal from the third and fourth source samples; 
a selected multiplier operatively coupled to the selected adder and configured to weight the sum signal 
with the selected weight to form a weighted sum signal; and 
is a third resulting sample builder operatively coupled to the selected multiplier and configured to produce the 

third resulting sample using the weighted sum signal. 

17. The signal processor of Claim 16 wherein the third source sample is the first source sample. 

20 18. The signal processor of Claim 16 wherein the fourth source sample is the second source sample. 

19. A computer system comprising: 

a computer processor; 
25 a memory operatively coupled to the computer processor; 

a signal processor stored in the memory and including at least one computer instructions which are executed 
within the computer processor to derive at least one resulting samples of a resulting signal from at least two 
source samples of a source signal, the signal processor including: 

30 at least five weights stored in the memory each of which is associated with one of the at least two source 

samples, 

wherein the respective weight associated with a first of the source samples in the derivation of a first 
of the at least two resulting samples is substantially equivalent to the respective weight associated 
35 with a second of the source samples in the derivation of a second of the at least two resulting samples, 

and 

further wherein the respective weight associated with the second source sample in the derivation of 
the first resulting sample is substantially equivalent to the respective weight associated with the first 
source sample in the derivation of the second resulting sample; 

40 

a composite weight builder operatively coupled to the at least five weights and configured to form first and 
second composite weights of the respective weights associated with the first and second source samples 
in the derivation of the first resulting sample; 

an sum adder configured to form a sum signal from the first and second source samples; 
45 a difference adder configured to form a difference signal from the first anc second source samples; 

a sum multiplier operatively coupled to the sum adder and configured to weight the sum signal by the first 
composite weight; and 

a difference multiplier operatively coupled to the difference adder and the composite weight builder and 
configured to weight the difference signal by the second composite weight. 

so 

20. The computer system of Claim 19 wherein the at least five weights include a selected weight associated with a 
selected one of the source samples in the derivation of a third of the resulting samples and another weight equiv- 
alent to the selected weight and associated with the selected source sample in the derivation of a fourth of the 
resulting samples; 

55 

further wherein the signal processor further comprises: 

a selected multiplier configured to weight the selected source sample with the selected weight to form a 
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selected weighted signal; 

a third resulting sample builder operatively coupled to the selected multiplier and configured to produce the 
third resulting sample using the selected weighted signal; and 

a fourth resulting sample builder operatively coupled to the selected multiplier and configured to produce 
the fourth resulting sample using the selected weighted signal. 

21 . The computer system of Claim 20 wherein the third resulting sample is the first resulting sample. 

22. The computer system of Claim 1 9 wherein the at least five weights include a selected weight associated with a third 
one of the source samples in the derivation of a third one of the resulting samples and another of the weights equiv- 
alent to the selected weight and associated with a fourth one of the source samples in the derivation the third result- 
ing sample; 

further wherein the signal processor further comprises: 

a selected adder configured to form a sum signal from the third and fourth source samples; 

a selected multiplier operatively coupled to the selected adder and configured to weight the sum signal 

with the selected weight to form a weighted sum signal; and 

a third resulting sample builder operatively coupled to the selected multiplier and configured to produce the 
third resulting sample using the weighted sum signal. 

23. The computer system of Claim 22 wherein the third source sample is the first source sample. 

24. The computer system of Claim 22 wherein the fourth source sample is the second source sample. 

25. A distribution system for distributing code (i) which is stored on a computer-readable medium, (ii) which is execut- 
able by a computer, and (iii) which includes at least one module, each of which in turn is configured to carry out at 
least one function to be executed by the computer, the distribution system comprising: 

at least five weights each of which are associated with one of at least two source samples, 

wherein the respective weight associated with a first of the source samples in the derivation of a first of at 
least two resulting samples is substantially equivalent to the respective weight associated with a second 
of the source samples in the derivation of a second of the at least two resulting samples, and 
further wherein the respective weight associated with the second source sample in the derivation of the 
first resulting sample is substantially equivalent to the respective weight associated with the first source 
sample in the derivation of the second resulting sample; 

a composite weight builder operatively coupled to the at least five weights and configured to form first and sec- 
ond composite weights of the respective weights associated with the first and second source samples in the 
derivation of the first resulting sample; 

an sum adder configured to form a sum signal from the first and second source samples; 
a difference adder configured to form a difference signal from the first anc second source samples; 
a sum multiplier operatively coupled to the sum adder and configured to weight the sum signal by the first com- 
posite weight; and 

a difference multiplier operatively coupled to the difference adder and the composite weight builder and config- 
ured to weight the difference signal by the second composite weight. 

26. The distribution system of Claim 25 wherein the at least five weights include a selected weight associated with a 
selected one of the source samples in the derivation of a third of the resulting samples and another weight equiv- 
alent to the selected weight and associated with the selected source sample in the derivation of a fourth of the 
resulting samples; 

further wherein the distribution system further comprises: 

a selected multiplier configured to weight the selected source sample with the selected weight to form a 
selected weighted signal; 

a third resulting sample builder operatively coupled to the selected multiplier and configured to produce the 
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third resulting sample using the selected weighted signal; and 

a fourth resulting sample builder operatively coupled to the selected multiplier and configured to produce 
the fourth resulting sample using the selected weighted signal. 

5 27. The distribution system of Claim 26 wherein the third resulting sample is the first resulting sample. 

28. The distribution system of Claim 25 wherein the at least five weights include a selected weight associated with a 
third one of the source samples in the derivation of a third one of the resulting samples and another of the weights 
equivalent to the selected weight and associated with a fourth one of the source samples in the derivation the third 
10 resulting sample; 

further wherein the distribution system further comprises: 

a selected adder configured to form a sum signal from the third and fourth source samples; 
is a selected multiplier operatively coupled to the selected adder and configured to weight the sum signal 

with the selected weight to form a weighted sum signal; and 

a third resulting sample builder operatively coupled to the selected multiplier and configured to produce the 
third resulting sample using the weighted sum signal. 

20 29. The distribution system of Claim 28 wherein the third source sample is the first source sample. 

30. The distribution system of Claim 28 wherein the fourth source sample is the second source sample. 

31 . The method of Claim 1 wherein the source signal is a recorded seismic signal. 

25 

32. The method of Claim 1 wherein the source signal is a video signal. 

33. The method of Claim 1 wherein the source signal is an audio signal. 

30 34. The signal processor of Claim 13 wherein the source signal is a recorded seismic signal. 

35. The signal processor of Claim 13 wherein the source signal is a video signal. 

36. The signal processor of Claim 13 wherein the source signal is an audio signal. 

35 
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